The effects of Zn Impurity on the Properties of Doped Cuprates in the Normal State 
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We study the interplay of quantum impurity, and collective 
spinon and holon dynamics in Zn doped high-T c cuprates in 
the normal state. The two-dimensional t-t'-J models with one 
and a small amount of Zn impurity are investigated within a 
numerical method based on the double-time Green function 
theory. We study the inhomogeneities of holon density and 
antiferromagnetic correlation background in cases with differ- 
ent Zn concentrations, and obtain that doped holes tend to 
assemble around the Zn impurity with their mobility being 
reduced. Therefore a bound state of holon is formed around 
the nonmagnetic Zn impurity with the effect helping Zn to 
introduce local antiferromagnetism around itself. The incom- 
mensurate peaks we obtained in the spin structure factor in- 
dicate that Zn impurities have effects on mixing the q=(7T, 7r) 
and q=0 components in spin excitations. 
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The effects of divalent transition metal Zn substitute 
for Cu in CuC>2 plane present much valuable informa- 
tion in understanding the mechanism of high tempera- 
ture superconductors. Zn 2+ has a closed d shell with spin 
s = 0, and acts as a very strong scattering center. As a 
result, the spin configurations and the electronic struc- 
tures around the nonmagnetic impurity Zn are strongly 
disrupted in both the normal state and superconduct- 
ing state. Below T c , there have been many experimen- 
tal [1-5] and theoretical [6-10] investigations to discover 
how the d-wave superconductivity is destroyed and what 
is the microscopic mechanism behind. As the normal 
state properties are more fundamental, some experimen- 
tal [11-15] and theoretical [16-18] works have been done 
to study the effects of nonmagnetic impurity Zn on the 
properties of normal state. It is believed that these stud- 
ies are of great help for understanding the peculiar be- 
haviors of normal state. Moreover, they shed light on 
the recent striking issue of the normal state pseudogap 
[19,20]. So far, there has been no clear picture of how 
the impurity interplays with the strong correlation back- 
ground. It still remains unanswered why the Zn impurity 
produced very strong scatter and how the impurity band 
forms in real space. 

In this paper, we perform numerical calculation to 
study the effects of Zn impurity on its surrounded Cu 
ions in the normal state. We want to find out how the 
Zn impurity influences the hole distribution and AF cor- 
relation background, which may bring about a better 
understanding of the fundamental relation between spin 



and hole. We start from the two-dimensional (2D) t-t'- 
J model and use fermion-spin theory [21]. Fermion-spin 
theory is based on the charge-spin separation, in which 
the single occupied constrain of t-t'-J model could be 
treated properly even in the mean field approximation. 
Within an improved Green function theory [22], we per- 
form numerical calculation for cases with only one Zn 
impurity and a small amount of Zn impurity, and the ef- 
fect of Zn concentration on some properties of the normal 
state are discussed. 

The essential physics of high-T c cuprates is well de- 
scribed by the t-J model on a square lattice. In the condi- 
tion that Zn substitutes Cu in Cu02 plane, we can model 
Zn impurity as vacant site, which has no coupling with 
surrounded Cu sites. We add the next-nearest-neighbor 
hopping term in our model to reproduce the realistic 
band structure, and start our study from the following 
Hamiltonian 

H = -t J2 (ClC ja + h.c.)-t' ]T {ClC Va + h.c.) 
-A^C&CW + J £ Si-Sj, (1) 

where (i, j) and (i, i') mean the summations over nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) pairs 
respectively, and I represents the site occupied by Zn im- 
purity. In our model the direct hopping among a Zn 
impurity and its surrounded Cu sites is forbidden. In ad- 
dition, to eliminate the doubly occupied sites of Cu ion, 
we introduce constraint C\ a Ci a < 1 for each Cu site. 
Also we introduce J2 a C\ Ci„ — 2 for Zn 2+ ion since it 
has a closed d shell. Therefore, the total number of elec- 
trons satisfies J^ia CiaCia = 2V — 2V/, + Nzn with 2V/, and 
Nzn representing the number of hole and Zn impurity, 
respectively. As the strong electron correlation manifests 
itself by the local constraint, the key issue is how to treat 
the constraint properly. 

Here we study the t-t'-J model within the fermion-spin 
theory [21] based on the charge-spin separation. We in- 
troduce Cjf = h\S~ and = h\Sf , where the spinless 
fermion operator hi describes the charge (holon) degrees 
of freedom, while the pseudospin operator Si describes 
the spin (spinon) degrees of freedom. Thus the low en- 
ergy behavior of the t-t'-J model (1) can be written as 

H = -t M + hjht)(SfS. + S7S+) 
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Fi and F 2 are the second- 
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where Jfj = [(1 - n^)(l - rij) - 4>\j]J with n\ repre- 
senting the hole concentration at site i and faj — (h\hj) 
being the order parameter of holon. Here we also in- 
troduce spinon correlation functions Xij = (S^ Sj~) and 
Xij — (SfSj). 

We introduce three double-time Green functions 

G h (i-j, T -T>) = -i6(T-T>)([h i (T); h+(r')}) 
= {{h i {T);hf{T'))) 
T-T>) = - i 6(T-T>){[S+: S-(r')]) 



Dl{i-j, T-T') = -ie{T-T'){[SU S z Ar')]) 
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= ((S z (t);S z (t'))), 

where Gh describes the behaviors of holon, and D s and 
D z describe the behaviors of spinon. Since the lattice 
translational invariance is not presented in cases with Zn 
impurities, we evaluate the equations of motion of the 
above Green's functions in real space. The double-time 
Green function ((A; B)) satisfies 

cj((A; B))„ = ([A, B] T ) U + (([A, H]; B)) u , (4) 

thus we could obtain the equations of motion of Gh 

(w - n)G h (i - j) u -2t^2 Xi,i+v G h{i + V~ .7% 
n 

-2t'Y J Xi,i+rG h {i + r-j) UJ = 5{i-j). (5) 

r 

We introduce Gh , a TV 2 x TV 2 elements matrix, to express 
the holon Green functions for a square lattice with N x N 
sites. And we could rewrite Eq. (5) as 



(lu - fi)G h -hG h = I, 
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where matrix h is decided by the NN and NNN spinon 
correlation functions, and I is an identity matrix. 

Based on Eq. (3) and (4), we also obtain the equations 
of motion of spinon Green functions D s and D z 

(jD s (i-j) w = 2^2Jl{l r) {e iti+v F 1 (i,i + r); j) u 
v 

~F\ (i + rj, i; j%} 
+8t'^2<t> ii i +T F 1 (i,i + t; j) u 

T 

uD z {i - j) u = J2 JtUn^i+vW, * + V, j) u 
n 

-F 2 {i + r],i; j%} 
+4t , ^2<t> itt+T {F 2 (i,i + t; j) u 

r 

-F 2 (i + T ,i; j) u }, (7) 



order spinon Green functions which are defined as 

W; j) u = ((S?S+; Sr)) u 

F 2 (i,l; j) u = ((S+Sr; S z )) u . (8) 

Going a step further, we establish the equations of motion 
of the second-order spinon Green functions 

f>Fi(i, I; j) u = 2Xi,Al - ]) - Xi,iS(i - j) 

+ (({^2[ 2 Jlj+r,( e lJ+V S i S l S t+ V ~ S i S l+v S f) 

n 

T 

+ 2<fii,i+rSi Sf S^_ T ]}; St)) u 
u)F 2 {i, 1; j) u = Xt,i$( l - j) ~ Xi,l$(i - .?) 

+ (({y][2Jl.l+n(St &l &1+T) ~ e l,l+riSt S l+ri Si) 
V 

+ 8t'J2[^,i + rSt +T SfSt 

T 

-cp hl+T s?sf +T snh (9) 

To the third-order spinon Green functions in the right 
hand side of Eq. (9) , we perform the improved decoupling 
scheme as described in Ref. 22, for example 

«Sf Sf 5+ ; ; SJ)) -+ ai {SlSl) ai ((S+- SJ)). (10) 



Therefore, the second-order spinon Green functions F\ 
and F 2 can be expressed by the Green functions D s and 

D z as 

I; j) u - T° + T\D(i - j) + r*D(Z - j) 

+ ^{r 2 ^^ + n - j) + TlD(l + rj - j)} 

n 

+ ^{rfi>(t + t- j) + t\d(i + t- j)} 

T 

u>F 2 (i, I; j) u = T° + T\D z {i - j) + Y\D Z {1 - j) 

+ ^{r 2 ^(» + rj - j) + TlD z {l + v - j)}, 



(11) 



where 



T° 1 =2xl l 6(l-j)-Xi,i6(i-j) 

r i = J i,i+r, e i,i+ri a i+riXi+v,l a l 
V 

+ 4< ' X! ^i.i+r^+TXi+r^m 

T 

r 2 T e ff 

r i = -Ji,ii v e i,i+r, a iXi,iai 
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T'i = -4*' ' 4>i,i+ T aiXi,iai 
rf = 8t'4>i t i +r aixliai 

?2 = XiM ~ 3) ~ Xi,iS(i - j) 

+8t' fiiA+Tpi+rXi+TjPl 

T 

?l = -2Jl{l v /3 iXi ,i(3i 
*2 - 2J?Jl v l3 i Xi,i0i 

(12) 

and 

St' ^ <l>l,l+T0iXi,l+T0l+T- 
T 

We also introduce two iV 2 x N 2 elements matrices D s 
and Dg to express the spion Green functions. Based on 
Eq. (7) and (11), we could obtain that 

lo 2 JJ s ~JDs=Cs 
u, 2 D*-S* D*=C*, (13) 

where matrices S, S z , C s and C z s are decided by the 
density of holon rii, and the spin correlation functions. 

We establish the self-consistent equations based on Eq. 
(6) and (12) to determine the correlation functions of 
holon and spinon, and also the vertex correction pa- 
rameters. Under the periodic boundary conditions, we 
have performed numerical calculations for 16 x 16 and 
20 x 20 lattices with different Zn concentrations respec- 
tively. Based on the experimental results of BSCCO near 
optimal doping [1,9], the parameters of t-t'-J model are 
taken as t/J — 2.5 and t' /t = —0.4 in our calculations. 

Firstly, we study the 20 x 20 lattice with only one Zn 
impurity in the optimally doped region (Sh = 0.15). The 
spacial distribution of holon density is calculated and our 
numerical results are shown in figure 1. We find that the 
holon density closet to the Zn impurity oscillates strongly, 
and the fluctuation diminishes rapidly away from the Zn 
impurity. The maximum density is obtained at sites two 
lattice distances away from the Zn impurity, and is about 
15% higher than the minimum density. Our numerical re- 
sults suggest that doped holes form a local region around 
the Zn impurity, whose size is about eight lattice cells as 
shown in figure 1. We also obtain that the magnetic 
modification introduced by Zn is mainly in the vicinity 
of Zn impurity. The NN correlation functions near the 



isolated nonmagnetic impurity in the undoped case have 
been discussed carefully in Ref. 22. We have found that 
Zn impurity strongly modifies the spin excitations, espe- 
cially the magnetic properties of the neighbor Cu. As a 
result, the AF correlation functions at the bonds closed 
to the impurities are enhanced. In the optimally doped 
region, we also obtain that Zn impurity enhances the NN 
AF correlation functions of spinons close to it. More- 
over, The doped holes have the effect to strengthen the 
AF correlations near the Zn impurity. The 63 Cu NMR 
study of YBa2(Cuo.9gZn . 01)306. 7 also find that the AF 
correlations are enhanced, not destroyed, around Zn im- 
purities [13]. Since the quantum fluctuation of spinons 
close to the nonmagnetic impurity is reduced obviously, 
we can divided the system into strong AF correlation re- 
gion and weak AF correlation region. The tendency of 
doped holes to assemble around the Zn impurity could 
rationalize anomalous charge localization effect, and the 
mobility of those holons closed to the Zn impurity could 
also be reduced. Therefore, a bound state of holon [23] 
is formed around the Zn impurity. Bound state of im- 
purity in normal state has also been predicted by the 
self-consistent T matrix approach [18]. 

We also study the cases with several Zn impurities in 
the optimally doped regime, and find that the holes play 
different roles in the strong AF correlation region and the 
weak AF correlation region. Our numerical results show 
that the AF correlations of bonds far from the Zn impu- 
rity reduces remarkably as hole is added into the systems. 
On the contrary, the bound state of holon has the effect 
to enhance the AF correlations around Zn impurity, and 
helps Zn to introduce local antiferromagnetism around 
itself. As Zn concentration increases, the Cu02 plane 
becomes a inhomogeneous mixture of strong AF correla- 
tion regions and weak AF correlation regions. And the 
doped holes tend to assemble at the strong AF correla- 
tion regions. 

We show the density of state of holon in cases with 
several Zn impurities in figure 2. Jung et al [24] have 
examined some samples and prove that Zn is uniformly 
distributed. In some configurations when Zn is uniformly 
distributed, our calculation show that due to the strong 
coupling between impurity and the conduction band, the 
width of the holon band decreases as Zn concentration 
increases. We find that the inhomogeneity of the spinon 
background increases the density of state of holon near 
the fermi surface, and a resonance peak of impurity states 
is found to get broader and stronger as Zn concentration 
increases. Inelastic neutron scattering study for the opti- 
mal doped Lai. 85 Sr .i5Cui_j / Zn y O4 indicates that a new 
in-gap Zn impuriy state is introduced at low tempera- 
ture [25]. Nonmagnetic defect structures at the surface 
has also been found to create localized low-energy ex- 
citations in their immediate vicinity in Bi2Sr2CaCu20s 
by performing low-temperature tunneling spectroscopy 
measurements with a scanning tunneling microscope [4]. 
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Our calculations show that the impurity state can sur- 
vive above T c , which is in agreement with the theoretical 
prediction [18] as well. 

To study the effect of the Zn impurity on the spin 
background around it, we introduce the spin structure 
factor 



s i (k) = 2s?s?, 



k(i-j) 



(14) 
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Here i represents Cu site around Zn impurity. Zn impu- 
rity is a scatter which has a strong effect on the AF cor- 
relation background [26]. As Zn concentration increases, 
the CuOi plane becomes a inhomogeneous mixture of 
strong AF correlation regions and weak AF correlation 
regions. In figure 3, the numerical results of the NN Cu 
sites around Zn impurity are shown for 16 x 16 lattices 
with different Zn concentrations. We obtain that, in the 
pure case and Zn lightly doped case (Sz n < 0.01), the 
spin excitations are dominated by a magnetic resonance 
peak located at Qaf — (7r, 7r). As the Zn concentra- 
tion increases, this peak decreases and there appear two 
second-high incommensurate peaks as shown in figure 3, 
which results from the mixing of q = (tt, tt), q = (tt,0) 
and g = components in spin excitations introduced 
by the strong impurity scattering. In Bulut's study of 
susceptibility of Zn doped high-T c superconductors, the 
similar behaviors are also obtained [17]. In addition, we 
find that the distances of the incommensurate peaks from 
q = (tt, tt) increase with doping, and these peaks become 
broad and weak in amplitude with the increasing of Zn 
concentration. Meanwhile, as a result of the increasing of 
the disorder introduced by the Zn substitution on the Cu 
sites, the peak at q = (tt, tt) decreases gradually as the Zn 
concentration increases, and disappears when Sz n > 0.1. 
Thus the result is consistent with experimental results of 
Zn-doped high-T c cuprates [11]. 

In summary, we have studied the interplay between 
quantum impurities, and collective spinon and holon dy- 
namics in Zn-doped cuprate in the normal state. Within 
a numerical method based on the Green function the- 
ory, the inhomogeneities of holon density distribution 
and antiferromagnetic correlation background in two- 
dimensional t-t'-J model with Zn impurities are investi- 
gated. We obtain the real space shape of bound state of 
holon surrounding the nonmagnetic Zn impurity. We also 
find that the doped holes help Zn to introduce local an- 
tifcrromagnctism around itself. In the cases with a small 
amount of Zn impurities, the influence of Zn impurity on 
the antiferromagnetic correlation background is studied. 
The appearance of incommensurate peaks in spin struc- 
ture factor indicates that Zn impurity is a strong scatter 
center, which has an effect on mixing the q — (tt,tt), 
q = (tt,0) and q — components in spin excitations. 
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FIG. 3. The Spin structure factor of the nearest neighbor 
Cu sites around Zn impurity for 16 x 16 lattices with different 
Zn concentrations. 




FIG. 1 . The distribution of holon density in a 20 x 20 lattice 
with only one Zn impurity. 




co — fj. 

FIG. 2. The density of state of holon at 8z n — 0.0 (dotted 
line), 0.01(dashed line), and 0.04 (solid line). 
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